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Abstract. Networks generated from Voronoi tessellations of space are prototypes for topo- 
logically disordered systems. In order to assess the effects of random coordination on the 
statistical properties of energy spectra, we analyze tight-binding models for planar topolog- 
ically disordered systems. To this end, the networks are generated by a simple topological 
model covering a wide range of naturally observed random structures. We find that the 
energy-level-spacing distributions exhibit level repulsion, similar to the spacings in the ener- 
getically disordered Anderson model of localization in the metallic regime. 
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1 Introduction 

In most natural structures disorder and random coordination are the rule rather 
than an exception. Therefore the investigation of random tessellations has become 
a subject of growing interest Jl]]. In addition to the structures found in nature, the 
Voronoi construction, a generalization of the Wigner-Seitz construction, allows to 
generate tessellations from sets of arbitrarily distributed points |Q. The random 
Voronoi tessellation has served as the prototype of topologically disordered structures 
in physics j|. The aim of this work is a first study of the influence of random 
coordination on the properties of the electronic spectra in topologically disordered 
structures. Similar structures are expected in liquid metals or alloys. Here we analyze 
a tight-binding (TB) model with neighbor hopping on certain two-dimensional (2D) 
random graphs. 

2 Tight-binding model on a random graph 

A simple topological model introduced by Le Caer allows to generate efficiently 
random mosaics that are very similar to Voronoi tessellations generated from dis- 
ordered arrangements of repulsive particles Starting point of this model is the 
triangular structure represented in Fig. |l]a (dashed lines). Disorder is achieved by 
random flips of the diagonal bonds which are performed with probability p, yielding 
the most disordered structures for p = 0.5. As an example, the central bond has 
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Fig. 1: a): Ordered p = triangulation (dashed thin lines) and mosaic (solid bold lines). 
The numbers indicate the number of edges of a closed loop (mosaic), or, equivalently, the 
number of neighbors of each vertex (triangulation). b): The same structure after a single 
flip. 



been nipped in Fig. |l|b. The dual graphs, the mosaics, obtained by connection of 
neighboring triangles, are represented as solid lines in Fig. |l|. 
The TB Hamiltonian is given as 

N 

where i,j label the sites of the 2D random graph. The total number of sites is N, and 
c\ (ci) denotes a fermionic creation (annihilation) operator on site i and we always 
use periodic boundary conditions. The hopping matrix element between neighboring 
sites i and j is tij = 1 and zero otherwise. The eigenvalues and eigenstates of the TB 
Hamiltonian are calculated employing standard diagonalization methods and Lanczos' 
algorithm. 



3 Density of States 

For the regular structures with p = 0, the density of states (DOS) can be obtained 
analytically. The dispersion relation of the triangular structure reads 

ex(fe) = 2 [cos^) + cos(/c y ) + cos(fca,. + k y )] , (2) 

where k = (k x , k y ). In the infinite lattice, the DOS is then given as 

PT (E) = {^Y 1 J dk S(E - e T (fc)) (3) 

where the integration is performed over the first Brillouin zone, i.e., k x ,k y £ [— 7T, it). 
This integral can be expressed in terms of complete elliptic integrals of the first kind. 
The resulting DOS is shown in Fig. |^. Note the logarithmic singularity at E = —2. 
For the mosaic structure at p — 0, the dispersion relation is found to be 

e M (fe) - ±V3 + e T (fe). (4) 

The corresponding DOS 

p M (E) = \E\p T (E 2 -3) (5) 
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Fig. 2: Left: DOS for the ordered triangulation at p — (dashed line) and the disordered 
triangulation at p = 0.5 (solid line). Right: Same as in the left panel but for the mosaic. 



with logarithmic singularities at E = ±1 is also shown in Fig. It is symmetric 
about E — because the ordered mosaic is bipartite. 

Fig. H also shows the DOS of the most disordered structures with p — 0.5, obtained 
by averaging over 50 realizations on a system of = 60 x 60 sites. These DOS 
are rather smooth, apart from oscillations at large energies which also show up in 
finite-size data for the ordered structures. The logarithmic singularities are washed 
out by the disorder, and the negative-energy peak for the mosaic structure almost 
disappears. However, the positive-energy peak for the mosaic structure survives for 
p = 0.5. 



4 Energy level distributions 

We now turn our attention to the statistical properties of the eigenvalue spectrum 
of the Hamiltonian (|l|). Here, we focus on the level-spacing distribution (LSD) P(s) 
of the normalized energy spacings s. In energetically disordered systems the metal- 
insulator transition (MIT) in the 3D Anderson model of localization is accompanied 
by a transition of the LSD. In the metallic regime, one finds level repulsion with 
P(s) = -Pgoe(s) of the Gaussian orthogonal random matrix ensemble (GOE) ||, 
whereas in the insulating regime, there is level clustering described by Poisson's law 
Pp(s) = exp(— s). Level clustering is also expected for ordered structures. For states 
at the MIT, an intermediate LSD has been observed [Q. 

We have computed the LSD of the complete spectrum by properly unfolding the 
spectrum || and averaging over 50 realizations on a system of 60 x 60 sites. In Fig. 
|, we show the results for random graphs with p = 0.01 and p = 0.5 in comparison to 
-Pgoe(s) and Pp(s). Both the mosaic and the triangulation data show an enhanced 
level repulsion with increasing p. The topological structures for p = 0.01 are still 
rather close to the ordered system. Thus it is no surprise that the LSD in Fig. 3 
still shows reminiscences of level clustering. However, from an investigation of the 
system size dependence we expect that even the p — 0.01 data tend towards the 
-Pgoe(s) result for increasing N. We have also studied the behavior of other spectral 
correlations such as the spectral rigidity A3 Our results indicate that there are 
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still characteristic deviations from the GOE behavior at the presently considered 
system sizes, although the system size dependence again shows a tendency towards 
GOE. This is further corroborated by an analysis of the statistical properties of the 
eigenfunctions. Results will be published elsewhere. 

5 Looking at the eigenstates 

In Fig. |^ we show the probability density of a typical eigenfunction i/j(r) for the 
triangular structure at energy E « 1.5 close to the center of the band. In agreement 
with the results presented above, we find that the eigenfunction extends over the whole 
system. This state is of course not a Bloch state, but nevertheless may be viewed as 
a superposition of eigenstates of the ordered system at energies E' = E ± AE, where 
AE represents an energy level broadening due to the disorder. An eigenstate of the 
weakly disordered system in Fourier space should therefore exhibit a broadened Fermi 
surface when interpreting E as the Fermi energy fio|| . We also plot the probability 
density of the Fourier-transformed eigenstate V>(fe) in Fig. ^ and compare it to the 
sharp Fermi surface of the ordered structure. The latter can easily be computed 
from the dispersion relation (||). For the strongly disordered case at p = 0.5 we find 
an almost four-fold symmetry, whereas the ordered triangulated structure only has a 
reflection symmetry with respect to the line k x = —k y . This is to be expected because 
for p = 0.5 both orientations for the diagonals (cp. Fig. |l|) are equally probable. Thus 
the four-fold symmetry of the square lattice is restored on average. Similar results 
hold for the mosaic structures. 

6 Conclusions 

We have considered the spectral properties of a TB Hamiltonian defined on certain 
2D random graphs. Although the systems are not energetically disordered, the LSD 
exhibits level repulsion much as for the energetically disordered Anderson model of 
localization in the metallic regime || and also for TB models defined on quasiperiodic 
tilings However, more evidence, such as the spectral rigidity mentioned above, is 
needed before we can conclusively decide whether the level repulsion is indeed iden- 
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Fig. 4: Left: Probability amplitude of a typical eigenstate at E ~ 1.5 in the random 
triangulation with p = 0.5 and N = 100 x 100. Dark regions correspond to high probability 
amplitudes. Center: Averaged probability amplitude of Fourier-transformed eigenstates at 
p = 0.01 with E and N as in the left panel. The lines indicate the Fermi surfaces of the 
ordered structure. Right: Same as in the center but with p — 0.5. 

tical to the behavior of GOE. We remark that preliminary results for eigcnfunction 
statistics may indicate systematic deviations from the GOE results for the most dis- 
ordered mosaic structures. 

We thank F. Milde for programming help in the initial stages of this project. 
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